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ABSTRACT 


The present study attempts to explore the possibility of crack detection in beams, through multi- 
tone harmonic excitation and higher order FRF analysis. A transverse crack in a beam has been 
modeled as a bilinear single degree freedom system, exhibiting breathing characteristics. A 
breathing crack opens or closes depending on the direction of the vibration. The response of such 
a system is examined under single tone, two tone and three tone excitations. The difference 
between the excitation frequencies is kept such that their harmonics occur at the natural 
frequency of the system, so as to magnify small amounts of bilinearity, which could have 
developed in a system originally taken as linear. Development of an equivalent polynomial 
nonlinear form to represent bilinear stiffness, is suggested for better understanding of the 
response harmonics. A step-by-step procedure for crack detection is developed. Numerical 
simulation is carried out through MATLAB toolbox to illustrate the procedure. Experimental 
work carried out on a laboratory beam set-up is also reported. 
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CHAPTER 1 


INTRODUCTION 


Damage detection by means of non-destructive testing plays an important role in ensuring the 
integrity of machine elements and structures. Vibration testing is an effective means of detecting 
irack development in structures. Fatigue cracks often exist in structural members that are 
subjected to repeated loading. Various studies have been carried out on the dynamic response of 
fatigue cracks, in an attempt to find viable vibration testing procedures for health monitoring. 
The crack models used in these analyses fall largely into two categories:(l) open crack models 
and (2) opening and closing or breathing crack models. 

1.1 Open Crack Model 

A beam is made more compliant by an open crack. The consequent reductions in the natural 
frequencies of free vibrations and changes in mode shapes have motivated a number of 
researchers to seek simple means for identifying and characterizing cracks in structures. Cawley 
and Adams (1979) developed an experimental technique to identify the cracks from changes in 
natural frequencies. Gudmunson (1982) used a first order perturbation method to predict the 
change in resonance frequencies of a structure due to cracks, notches or other geometrical 
changes. Christides and Barr (1984) examined the special case of symmetric double sided crack. 
They used variational principles to formulate a fourth order partial differential equation for a 
Bemoulli-Euler beam with a double sided crack. Mathematically and experimentally, they found 
that the drop in the first natural frequency was very small unless the crack was fairly large. 
Gounaris and Dinarogonas (1988) modelled a crack with an equivalent rotational spring 
constant, calculated via an application of Castigliano’s theorem. They also found that the 
changes in the natural frequencies were significant only for large cracks. Rizos et al (1989) 
conducted experiments to detect crack depth and location from changes in the mode shapes of 
cantilever beams. They forced a cracked beam with known crack characteristic in resonance and 
measured vibration amplitude at two arbitrarily chosen position on the beam. Further research on 
crack identification via natural frequency changes was done by Anifantis et al. They found that 
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he size and position of large cracks can be found if the change in the first two natural 
'requencies were known for a simply-supported Bernoulli-Euler beam. Rajab and Al-Sabeeh 
1991) found that knowledge of the changes in the first three natural frequencies of a 
rimoshanko beam was sufficient to determine the size and locations of large cracks. Gounaris et 
al (1991) modelled cracked hollow beams with the finite element method , but encountered a 
sensitivity problem with small cracks. Ostachowicz and Krawczuk (1991) examined the effect of 
more than one crack on the natural frequencies of a cantilever beam, yielding comparable results. 
Reference can also be made to the studies by Pandey et al (1991), Nandwana and Maiti (1997), 
Tsai and Wang (1996),on analysis of open crack models. 

1.2 Breathing Crack Model 

Open crack models, however are inadequate in correctly representing the phenomenon due to the 
assumption that the cracks are always open in vibration. A breathing crack model is a more 
realistic representation. The basic concept is that the crack alternatively opens and closes 
depending on the direction of vibration. When the crack is closed, the beam acts approximately, 
as a homogeneous beam with no crack, while when the crack is open, a local reduction of 
flexural rigidity occurs. Recently, increasing efforts have focussed on vibration analysis using 
opening and closing models to simulate a fatigue crack. A single-degree-of-freedom, non-linear 
model was proposed by Friswell and Penny (1992). A mathematical model of a beam, with a 
closing crack was developed by Krawczuk and Ostachowicz (1994), who examined the 
possibility of crack detection on the basis of changes of the high harmonics in the frequency 
spectrum. Abraham and Brandon (1995) presented a substructure model to predict the vibration 
properties of a beam with a breathing transverse crack. Dimarogonas and Papadopoulos (1983) 
have obtained analytical solutions for the closing crack, for rotating shafts, under the assumption 
of large static deflection. Gudmunson (1983) pointed out that the alteration in natural frequencies 
due to a real fatigue crack is much lower than the drop caused by the narrow notch. Ibrahim et al 
(1987) confirmed the result and employed a bond-graph model to simulate the dynamic behavior 
of a cantilever beam including a non-linear fatigue crack. Qian et al (1990) observed that the 
difference of the amplitude response between the integral and cracked beam is reduced if a 
closing crack model is considered. The salient features of the studies cited above can be 
summarized as - 
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) Natural frequency reduction for a fatigue crack (breathing crack) is much smaller than for 
an open crack and that fatigue cracks would be difficult to recognize by frequency 
monitoring. Also, crack detection by an open crack model would underestimate the crack 
severity if the crack was actually growing under fatigue loading conditions. 

ii) The experimental studies reveal the presence of super/sub-harmonic vibration 
phenomena. 

iii) Detection of fatigue cracks should be more reliably based on non-linear features of FRF, 
rather than the natural frequency shift. 

1.3 Present Study 

The focus of the present study is to analyze the bilinear model of a cracked beam. The bilinear 
model is discussed in the next chapter. The proposed crack detection procedure involves multi- 
tone harmonic excitation. Numerical simulation of the response is discussed in Chapter 3. An 
equivalent polynomial form of nonlinear modeling is done and the response is represented as a 
Volterra series, in order to rationalize the numerical results. Experimental investigations are 
described in Chapter 4. Chapter 5 gives the conclusions and scope for further investigations. 
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CHAPTER 2 


NONLINEAR MODELS OF A CRACKED BEAM 

Nonlinear models of a cracked beam are discussed in this chapter. Damping is assumed to be 
linear and only stiffness nonlinearity is considered. Two forms of nonlinear spring force 
representation are considered - (i) Bilinear function and (ii) Polynomial function. 

Bilinear oscillators, due to their simplicity, appear to be appropriate to simulate the behaviour of 
fatigue cracks in structures. The fundamental assumption is that the beam vibrates mainly in its 
first mode. The exactness of this assumption depends on the position and the kind of force 
exciting the structure. A transverse crack, as shown in Fig. 2.1, will alternately open and close 
under the action of an excitation force. When the crack is in the compression (closed) the 
equivalent stiffness may be regarded as that of a beam without a crack, due to the contact of the 
crack walls. On the other hand, when the crack is in tension (crack open) the beam stiffness will 
reduce near the crack. The equivalent stiffness, of this model, therefore takes two different 
values depending on whether or not the crack is open. The equivalent damping of the beam is 
assumed to be constant although the system has increasing damping due to the friction forces at 
the point of the crack closure. Such a model was considered by Shen and Chu (1992), by 
deriving a bilinear equation of motion for each vibration mode of a simply supported beam. 
However, an analytical representation of the response of a bilinear system is approximate. 
Galerkin procedure was employed by Shen and Chu to simulate the response. Rivola and White 
(1998) idealized the cracked beam as a simple SDOF system with bilinear stiffness 
characteristics. They presented a procedure to obtain approximate analytical expressions of the 
response of bilinear oscillators. This procedure is briefly outlined below. 

2.1 The Bilinear Model 

The simple bilinear single-degree-of-freedom bilinear model of a cracked beam, shown in 
Fig.2.2, can be described by (Rivola and White, 1998) 

mil + cii + k (u)u = f (t) (2-1) 
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where u(t ) is the displacement of the oscillator, m is the mass, c the damping coefficient and 
k («) is a piecewise stiffness function defined as 
k ( u ) = {k,if u < 0 

kaif u > 0 } with 0 < a < 1 (2.2) 

The relationship between the restoring force and the displacement of the bilinear system is 
depicted in Fig.2.3. For free undamped vibration of the oscillator, the period of the system is 
composed of two half sinewaves, and a bilinear radian frequency co 0 , is defined as (Choi and 
Noah, 1988) 

co 0 = 2 co x co 2 l(co x + a> 2 ) (2.3) 

where 

a> ] = 4kTm co 2 = yfak/m (2.4) 

Obtaining a closed form solution of equation (2.4) is quite difficult. Kraczuk and Ostachowicz 
(1994) found the solution for a harmonic excitation force, by approximating the stiffness force 
through a square-wave function s(t), having a frequency equal to the forcing frequency as 
follows 

mil + cu+ k[ 1 - (1 - a)s(/)]« = / ( t ) (2.5) 

The square- wave function s(t), given in Fig. 2.4, changes between 0 and 1, and can be 
approximated by the complex Fourier series 

(2 - 6) 

/J *i—0O 

where co is the excitation frequency and q„ are Fourier coefficients 

q n =^— ^cos[— (« + l)] if n*0, q n = 1/2 ifn = 0 (2.7) 

nx 2 

Using equations (2.6) and (2.7) in (2.5), the equation of motion becomes 

mil + cu + k[ 1 - (1 — a)^q„e lnoM ]u = e~ ia * + (2.8) ,| 

-oo 2 2 

It is to be noted that the above equation represents only an approximation of the spring force 
variation. Rivola and White (1998) further employed the harmonic balance method to transform 
the above equation into a set of linear equations, which were then solved numerically. 
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2.2 The Polynomial Model 


While, representation of a cracked beam through a bilinear spring model helps in understanding 
the role of the opening and closing of crack, on the spring forces developed by the beam over a 
cycle of vibration, the absence of an analytical solution restricts its utility in understanding the 
response phenomenon. In this context, an interesting study reported by Sundermeyer and Weaver 
(1995) can be discussed. They carried out a numerical simulation on a bilinear oscillator. It was 
excited by two sinusoidal forces with different forcing frequencies. They observed that if the 
difference of the two forcing frequencies is approximately equal to the fundamental frequency of 
the system, in addition to the peaks at the two forcing frequencies, a peak is also observed at the 
natural frequency, in the FRF of the response of the cracked beam. This has been graphically 
described in Fig.2.5. Sundermeyer and Weaver, while pointing at the possibility of employing, 
the occurrence of this additional peak at the natural frequency, as a tool to detect the presence of 
a crack, do not provide analytical reasons for this behaviour of a beam with crack. 

Analysis of the above behaviour can be attempted through a polynomial representation of the 
cracked beam stiffness. The emphasis of the present study is to investigate the response pattern 
of the bilinear oscillator through its better analytical representation, which can help in modifying 
and refining the crack detection procedure suggested by Sundermeyer and Weaver. 

The nonlinear spring force described by the bilinear oscillator, representing a cracked beam can 

be approximated through a polynomial function as illustrated in Fig. 2.6 and given by 

g(x) = k x x + k 2 x 2 + k 3 x 3 + V 4 + M* C 2 - 9 ) 

This polynomial function approximation renders the governing equation (2.1 and 2.2), as 
mx + cx + k x x + k 2 x 2 +k 3 x 3 +k A x 4 +k 5 x 5 = fit) (2.10) 

Using Volterra series representation (Schetzen, 1980) of the response, x(t) is expressed as 

*(0 = 2 X (0 

n»l 


( 2 . 11 ) 



( 2 . 12 ) 


with the nth order response component being given by 

oo oo oo 

X n (0 = f J • fh„(r l ,T 2 ,...,T„)f(t-T,)f(t-T 2 )...f(t-T„)dT 1 dT 2 ..JT„ 


h n (ri,r 2 ,—,T n )is nth order Volterra kernel and its Fourier transform provides the nth order 
Frequency Response Function (Schetzen, 1 980) as 


H n (£D„o) 2 ,...,a> n )= | \ h n ( r„r 2 ,.. 


;T n )fle- Ja,,Tl dT l dT 2 ...dT n 

/=! 


(2.13) 


Employing the above, the individual response components (equation 2.12), can be expressed, 
after some algebra as 

x n (t) = (A/2 )” E n C q Hr{co)e j “ 0 <p<n; 0 <q<n (2.14) 

p+q=n 


where the following brief notations have been used 

= ,-cu) co pq =(p-q)a> (2.15) 

V "" “ V ^ V "" ‘ mm v J 

p times q times 


2.2.1 Response to Single-Tone-Excitation 

Restricting the spring force nonlinearity, in equation (2.10), upto the cubic term only and a single 
tone harmonic excitation 


f(t) = A cos cot = — e jM + — e~ jM 
2 2 

the Volterra series response becomes 
x(t) = x, (0 + x 2 (0 + x 3 (0 + .... 

where the individual response components *,(/), * 2 (0> x 3 (0 ^ 


x \(t)= 

~a 

4 


— e> a '( <_r ) +—e 


dt 


-e jM hAr)e- jM +-e- Ja *h l (r)e Jal 
2 2 




■e J * a Jh l ;(T)e- J *dr + -e- J * jh,(r)e JM dr 


= - e ,M H x (a) + -e~ jM //, (-<y) 
<2 2 

= |//,(^)e^+|/f,(-cy)e-^ 


(2.16) 

(2.17) 


(2.18) 
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* 2 (0= J J M r i.r 2 ) 


^_ e jo>(l-T X ) + ^_ e -j»>{l-r i) 
2 2 


ii g ) + A e -JaV-Tz ) 

2 2 


dx x dx 2 


a oo 


= | | /i 2 (r,,r 2 ) 

-a 

A 


- a -oo 
#2 


> ./w(2/-~r, -r 2 ) 


, e - jo ) r ^ jmr 2 _ e J '» T l - jt » r 2 j A c ~ lj ( X + jm ^ ]< 0T r Ly 




^( 2 1, 


4 —oo -a 4- -a 

+ ; 7 l|/> 2 ((,,r 2 )e y "''-'" r ‘*,dr a +^-e~ 1J “ 5jfc 2 «„r 2 )e'“'***'>rfr 
4 -a 4 -a 

= ^V 2< “// 2 (<»,<») + ^- e - j °*H 2 (a -a) + ~e J 2 °*H 2 (-co,co) + ~e~ jat H 2 (-co-co) 

A 2 . 

= — e j2a * H 2 (co, c o ) + — H 2 (co-co) + conjugate terms 
4 4 


2.19 


* 3 (0= I JI^ 3 (r„r 2 ,T 3 ) 


—e J ‘» (,-r *) + —e~ J<a( -‘~ T ^ 

* 

^ g J«('-r 2 ) + j£ g 

_2 2 


.2 2 




.2 

= iPsO: 


j. T 2 » T 3>[ |- 


dr l dx 2 dz 3 


. e jM 2 r-Ti-r 2 ) Q~j o>r \ + j° >t z | _ e J m \-i<» T z 


a 2 

+ g-ljax+jtOTx+jmi 

4 


j£ g M'-rj ) + _f£ g -J<B(<-r 3 ) 


4 

dz x dx 2 dx 3 


= fJjA3(r l5 r 2 ,r 3 )[^e >(3 '- r, - r2 - r3) +^e >( '- Tl - r2+r3) + £_*><*-*»+**-«*> 


[ 21 g ;<B(-/-ri+r 2 +r 3 ) | ^ | ^ c ->ai(-f+r|-r 2 +r 3 ) 

8 8 8 

+ £ e M-<«, + r J -, J ) + £ e >(- 3l «, t , J+ r ) )j (/ , . 

8 8 j 1 2 3 

A 3 0 A 3 

= — e j3a * H 3 ( a ), co, c 0 ) + e Ja * H 3 (co, co-co ) 

8 8 

■ZA 3 A 3 

+ — e~ jM H 3 (co -co-co) + — e~ j3a * H 3 (-co-co ,-co) 

8 8 

A 3 3 A 3 

= — e j3a * H 3 (co,co,co) + e JC *H 3 (co,co ,-co) + conjugate terms 

8 8 


( 2 . 20 ) 
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In equations (2.18)-(2.20), the term H x (a> x ) is the first order Frequency Response Function of 
the system and can be readily shown to be given by (Schetzen, 1980) 

= . ( 2 . 21 ) 

‘\j(k x -mco 2 ) 2 + (ccy)~ 


and the higher order FRFs, H 2 {co x ,co 2 ) and H 3 (co x ,co 2 ,a> 3 ) can be expressed in terms of the first 
order FRF H x {a> x ) as given below (Chatterjee and Vyas, 1999) 

H 2 (®, , 0 ) 2 ) = —k 2 H x (co x )H X ( co 2 )H X (co x + co 2 ) 


Hy (©,, co 2 , co 3 ) = H x (co x )H x {co 2 )H x (a> 3 )H x (o) x +co 2 + co 3 )* 



{H x (o) x +q) 2 ) + H x (co 2 +(Oy) + H x (co 3 +co x )}-k z 


( 2 . 22 ) 


Collecting terms of identical frequency, in x x ,x 2 ,x 2 from equations (2.18) - (2.20), gives the 
overall response x(t ) in terms of its various frequency components as 

x(t) = z 0 + 1 z(<y, )| cos(ry, t + (f > x ) + \z(2co x )| cos(2® x t + (j ) 2 ) + \z(3a > x )| cos(3®, / + ^ 3 ) + ... (2.23) 

where 


3 

z(co x ) = AH X (o,) + — A 3 Hy(a> x , a> x ,-o) x ) 

+ — A 5 H 5 (G) x ,a> x ,a) x -(v x ,-co x ) + higher order terms 
16 

A 2 A 4 

z(2a>, ) = — H 2 (a)„(V 1 ) + — H 4 (<v x ,a> l ,(V x -a >, ) + 

A 2 5 A i 

z(3ty, ) = — H 3 (<y, , co x , <y, ) + — — H s (o>, ,a) x ,a) x -co x -a , ) + . 
4 16 

3 

z 0 = A 2 H 2 ( co x ,-<y, ) + - A 4 H 4 (<u, , co, ,-<y , , <o x ) + 


(2.24) 


Restricting, to the first to third FRFs alone (terms H x (co x ) to H 3 (cy, , a> x ,-a> x ) , i.e. excluding 
H 4 (o, , co x -a > x , o) x ), H 5 (<y, -ft>, ) and other higher order FRFs), in equations (2.24), 

it can be noted, from equations (2.21) and (2.22), that the maximum response at frequencies 
co x , 2co x , 3 co x will occur respectively when these frequencies are equal to the linear natural 

frequency (^k x 7m). f 
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2.2.2 Response to Two Tone Excitation 
For a two-tone harmonic excitation, 

f(t) = Acosw,(t) + Bcos 0 2 t (2.25) 

of the system described by equation (2.10), (with spring force including terms upto cubic only), 
the various frequency components of the response x(t) can be shown to be 


3 3 15 

z(co , ) = AH, (co,) + -A 3 H 3 (co, , co, ,- co , ) + - AB 2 H 3 (a, , co 2 ,-co 2 ) + — A 5 H s (©, , <u, , a, - 0, -co , ) 

4 2 16 

+ — AB 4 H 5 (a, ,co 2 ,co 2 - co 2 ,-co 2 ) + ^~ A 3 B 2 H ^(co,,®, ,-co, co 2 ,-co 2 ) 

16 8 

3 3 15 

z(co 2 ) = BH, (co 2 ) + -B 3 H 3 (a 2 ,0 2 ,-a 2 ) + -BA 2 H i (0„0 2 ,-co 2 ) + — B 5 H s (co 2 , co 2 , <o 2 - co 2 ,-co 2 ) 

4 2 16 


15 


4- — BA H s (co l , co x -co x - co , -co 2 ) + — j B A H 5 (co 2 , cy 2 -co 2 9 co x -co x ) 

16 8 

z( 2 &, ) = A- H 2 (co, , a , ) + H A (0, , 0, , a>, - 0 , ) + ~A J 2 B 2 H A (0 , , <y, , co 2 , e> 2 ) + ... 

jL JL <Z» 

g 2 ^2 ^ 

z(2<y 2 ) = — tf 2 (®p< ) + — H a (co 2 ,0 2 ,co 2 ,-<o 2 ) + - A 2 B 2 H, (©, ,co 2 ,. a> 2 ) + ... 

z( 3£0.) = — (# p rw, ) + ^ yl 5 + f A 3 B 2 H 5 (0 x ,0 x ,0„0 2 ,-0 2 ) + . 

4 16 4 


B 3 5 5 

z(3ty 2 ) = — H 3 (a> 2 , o 2 , co 2 ) + — B 5 H 5 (co 2 , co 2 , co 2 , a> 2 ,-co 2 ) + - B 3 A Z H 5 (co 2 , co 2 , <u 2 , < 0 , - 0 , ) + ... 

4 1 6 4 

3 3 

z(<y, + © 2 ) = ABH 2 ( 0 ,, 0 2 ) + -A 3 BH 4 (©,,©, , <a 2 ) + ~ AB 3 H A ( 0 , , co 2 ,-co 2 ,~ 0 2 ) + . .. 

Z z 

3 5 15 

z(2<y, + co 2 ) = - A 2 BH 2 ( 0 , , co, , <y 2 ) + - A 4 BH 5 (oo x , 0 , , ( 0 , ,-co , , ) + ~ A 2 B 3 H } ( 0 , , co, , co 2 , a> 2 ,-co 2 ) + .. . 

4 4 O 

3 5 15 

z(2co, - co 2 ) = — A 2 BH 2 ( 0 „ 0 „- 0 2 ) + -A a BH s ( 0„0,,0, ,-co, ,- 0 2 ) + -A 2 B 3 H,(0,,0„0 2 ,-0 2 ,-co 2 ) 4 

4 4 o 

3 5 15 

z(lco 2 + 0 ,) = -B 2 AH 2 ( 0 „ 0 „ 0 2 ) + -AB 4 H i ( 0 „ 0 2 , 0 1 , 0 2 ,- 0 2 ) + —A 3 B 2 H s ( 0 „ 0 „- 0 „ 0 2 , 0 2 ) + ... 

4 4 O 

3 5 15 

z(lco 2 - co , ) = — 5 2 AH 3 (-©, ,co 2 ,co 2 ) + -B 4 AH 5 (-co, ,co 2 ,co 2 ,co 2 ,-co 2 ) + — A 3 B 2 H t (co, ,-co, - 0 , , 0 2 , 0 2 \ 

4 4 o I 

z 0 =A 2 H 2 ( 0 „- 0 ,) + B 2 H 2 (0 2 ,- 0 2 ) + ^A 4 H a (0„0, ,- 0 , - 0 ,) 


D 4 22T /V,* xiv _/vy U 


(2.26) 



2.2.3 Response to Three Tone Excitation 
Similarly, for a three-tone harmonic excitation, 
f(t) = A cos co l (/) + B cos co 2 t + Ceos co 3 t 

the frequency components of the response can be worked out to be 


(2.27) 


3 3 3 

z(co x ) = AH X (co x ) + — A 3 H 3 (a > x , co x ,-co , ) + — AB 2 H 3 (a> x ,a) 2 ,-o) 2 ) +—AC 2 H 3 (po x ,a> 2 ,-a > 3 ) 

4 2 2 

+ ~ A ! H s (rU| , , £t>, -iy 1 ,-® I ) + -^ r ^S 4 // 5 (<a) l ,n; 2 ,<u 2 -<a 2 ,-£y 2 ) 

16 16 


+ — AC 4 if j (m, , , m 3 - co 3 ,-<u 3 ) + ^- ^4 3 5 2 ff 5 (< 2 >, , <y, , co 2 ,-<o 2 ) 

16 


15 
8 

1 ^ . 


+ — A 3 C 2 H s (co x , co x ,-co x , <y 3 ,-tu 3 ) + — AB 2 C 2 H S {co x , co 2 -a> 2 ,co 3 -co 3 ) + ... 

8 4 

3 3 3 

z{co 2 ) = BH X (a> 2 ) + -B 2 H 3 {co 2 ,a) 2 -co 2 ) + -BA 2 H 3 {co x -co x -co 2 ) + -BC 2 H 3 (g> 2 ,0) 3 ,-a> 3 ) 

+ — B 5 H 5 (co 2 ,co 2 ,co 2 -o) 2 ,-o} 2 ) + -^-BA a H 5 (co x ,(o x ,-co x -a> x ,a> 2 ) 

16 16 


+ — BC 4 H j (ru 2 , <a 3 , a> 3 - co 3 ,-co 3 ) + ~B 2 A 2 H 5 (co 2 ,Q) 2 ,-o) 2 , a> x ,- 0 ) x ) 


_5_ 

16 

15 

8 


8 


+ ^ B 2 C 2 {a 2 , co 2 -o) 2 , (o 3 ,-o) 3 ) + ~A 2 BC 2 H s (cd x -co x ,o) 2 ,co 3 -o ) 3 ) + . 


3 3 

z{a > 3 ) = CH X (<y 3 ) + - C 2 H 3 (<u 3 , a 3 -a > 3 ) + — CA 2 H 3 (co x -co x ,-oj 3 ) 

3 

+ — CB 2 H 3 (qj 2 ,-o ) 2 , co 3 ) + ... 

A 2 A 2 3 , , 

z(2cu, ) = —H 2 (co x ,g) x ) +—H 4 (cd x ,q) x ,co x ,-<a x ) + - A 2 B 2 H 4 (eo x , co x , a > 2 ,-(o 2 ) 
2 2 2 

3 

+ — A 2 C 2 H 4 (o) x ,co x ,q) 3 -£o 3 ) + ... 

z(2o) 2 ) =—H 2 (q) 2 ,o) 2 ) + ^-H 4 (co 2 ,o} 2 ,o) 2 -o) 2 ) + ^A 2 B 2 H 4 (q} x ,-(v v co 2 ,q} 2 ) 
2 2 2 
3 

+ - B 2 C 2 H 4 (q) 2 , co 2 , o) 3 ,-o) 3 ) + ... 

z( 2 g > 3 ) = — H 2 (o) 3 ,o) 3 ) + ^-H 4 (o} 3 ,qj 3 ,o} 3 ,-g ) 3 ) + | A 2 C 2 H 4 ((v x -oj x , co 3 ,<v 3 ) 


(2.28) 


+ ^ B 2 C 2 H 4 (a> 2 ,-oj 2 , ct) 3 ,> <n 3 ) + . 
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A 3 5 

z( 3 ft), ) = — H 3 (ft) , , ft), , ft) , ) + — 4 5 H s (o, , ft) , , ft), ,-G), ,-ft), ) 

4 16 

+ — 4 3 B 2 H. (ft), , co, ,o).,o) 2 ,-co 2 ) + — A 3 C 2 H. (ft), , co. , ft), ,a>, -co z ) + . 
4 4 


4 


z(3ft) 2 ) = ^ H 3 (o 2 , co 2 , co 2 ) + ^ B 5 H s (o 2 , co 2 , co 2 , co 2 ,- 0 ) 2 ) 


+ — B A H s (co 2 ,co 2 ,o) 2 ,co ] ,-co ] ) + ... 


C 3 5 

z( 3®3 ) = — H 3 (ft) 3 , <y 3 , co 3 ) + — C 5 H 5 (ft) 3 , co 3 , co 3 , co 3 -co 3 ) 

4 16 

+ — C 3 A 1 H Aco 3 ,co 3 ,co l ,co x -co x ) + — C 3 B 2 H s (co 3 ,co,,co 3 ,co 2 ,-co 2 ) + ... 

4 4 

3 3 

z(ft), + ft) 2 ) = ABH 2 (ft), , ft) 2 ) + — 4 3 BH A (ft), , ft), ,-ft), , ft) 2 ) + — AB 3 H 4 (a), , ft ) 2 ,-fi ) 2 ,-ft ) 2 ) 

2 2 

+ 4BC 2 // 4 (ft), , ft ) 2 , ft ) 3 ,-fi ) 3 ) + ... 

3 3 

z(ft ) 2 + ft) 3 ) = BCH 2 (ft) 2 , ft) 3 ) + - B 3 CH A (ft) 2 , ft ) 2 ,-ft) 2 , ft) 3 ) + — AB 3 H a (a), , fi) 2 , fi ) 2 ,-ft) 2 ) 


+ 34 2 BCtf 4 (ft), ,-ft), , ft) 2 , ft> 3 ) + ... 

3 3 

z(ft), + fi) 3 ) = ACH 2 (ft), , fi) 3 ) + - A 3 CH a (ft), , ft), ,-ft), , ft) 3 ) + - AC 3 H a (ft), , ft) 3 , ft ) 3 ,-ft) 3 ) 

2 2 

+ 34B 2 C // 4 (ft), , ft) 3 , o 2 ,-fi) 2 ) + .. . 

3 3 

z(ft), - fi ) 2 ) = 4B/7 2 (ft), ,-ft) 2 ) + - A 3 BH a (ft), , ft), ,-ft), ,-ft ) 2 ) + -AB 3 H a (co 1 , co 2 ,-ft ) 2 ,-fi ) 2 ) 

2 2 

+ 3 ABC 2 H a (ft), ,-fi) 2 , ft ) 3 ,-fi) 3 ) + ... 

3 15 

z(ft), + ft) 2 +ft) 3 ) = — ABCH 3 (e), , o 2 , ft) 3 ) + — A 3 BCH j (ft), , ft), ,-e), , ft) 2 , a) 3 ) 

+ — 4 B 3 CH j (ft), , ft) 2 , ft ) 2 ,-ft) 2 , ft) 3 ) + ^ ABC 3 H } (ft), , ft) 2 , ft) 3 , fi ) 3 ,-ft) 3 ) 

4 4 

+ higher order terms 


3 

z(2fi), + ft ) 2 ) = - 4 2 Btf 3 (ft), , ft), , fi ) 2 ) + 4 4 BH } (ft), , ft), , ft), ,-ft), , ft) 2 ) 

15 <2d3i 


5 44; 
4 


+ — 4 B //j (a), , ft), , ft) 2 , ft ) 2 ,-ft) 2 ) + - 

8 


z(2co ] -cd 2 ) = — A 2 BH y ( cd x , 6) x 9 -g > 2 ) + - AB 4 5 (o ) , , <y 2 , 2 r ) 

4 4 

+ — 4 2 B 3 # 5 (ft), , fi), , ft) 2 ,-ft ) 2 ,-ft) 2 ) + -^- 4 2 BC 2 H 3 (ft), , ft), , ft ) 2 , ft ) 3 ,-ft) 3 ) + . .. 

8 8 
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3 5 

z(2g ) 2 + co , ) = — B 2 AH 3 (a), , co 2 , a > 2 ) + — AB* H s (©, , co 2 , co , , a> 2 ,-co 2 ) 

4 4 

+ V A 3 B2H S » ®1 > _£a i , > ^2 ) + V 2C 2// S > <*>2 > <°2 > -®3 ) + - 

O O 

3 5 

z(2<a 2 - <», ) = - S 2 ^j (-co, , a) 2 , a) 2 ) + - B 4 AH 5 (-ft), , ft) 2 , fi) 2 , co 2 -co 2 ) 

+ T A 5 5 2 ^5 -^ 1 ,^ 2J ^ 2 ) + ^^5 2 C 2 // J (-^ 1 ,ffl 2 ,<y 2 ,<»3 ,-ft) 2 ) + ... 

O O 

3 5 

z( 2 fi > 2 + a) 3 ) = —A 2 CH i (a> i ,a) 2 ,co 2 ) + — A A CH s (ct) ] , co^co. , co., -&>,) 

4 4 

+ — A 2 B 2 CH 5 (a ), , co , , a ) 2 -co 2 , ft) 3 ) + — yt 2 C 3 H s (co, , co , , a) 3 , co 3 ,-co 3 ) + .. . . 

4 4 

3 5 

z(2co , -co 3 ) = — A 2 CH 3 (co |,o, ,-co 3 ) + — A 4 Ci/ 5 (a), , a), , ft), ,-a), ,-a) 3 ) 

4 4 

+ — A 2 B 2 CH 5 (ft), , £ 0 , , co 2 ,-co 2 ,-co 3 ) + — A 2 C 3 H s (a), , a), , co 3 , ft ) 3 ,-a) 3 ) + .... 

4 4 

z 0 =A 2 H 2 (co, ,-ft ), ) + B 2 H 2 (a ) 2 ,-a) 2 ) + C 2 /7 2 (ft ) 3 ,-ft) 3 ) 

3 3 

+ - A 4 // 4 (ft), , a;, ,-£ 0 , ,-ft), ) + 6 A 2 C 2 H 4 (ft), ,-ft), , a ) 3 ,-ft ) 3 ) + -B 4 H A (oo 2 ,co 2 ,-co 2 ,-o 2 ) 

4 4 

3 

+ 6B 2 C 2 H A (co 2 ,— ft) 22 , ft > 3 ,-ft) 3 ) + - C 4 H 4 (ft) 3 , a), ,-fi > 3 ,-ft) 3 ) + ... 

It can be seen, from equations (2.24), (2.26), (2.28), that for a system with polynomial type of 
nonlinearity (equation 2 . 10 ), excited by harmonic forces, the response comprises of components 
at frequencies, which are multiples and various linear combinations of the excitation frequencies. 
It is to be noted that, if the system is linear, the response would contain frequency components at 
the excitation frequencies £y,, a> 2 , co 3 only. This fact is employed in the proposed crack 

detection procedure, which is discussed in the next chapter. 
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Fig. 2.3 Piecewise restoring force 



Fig. 2.4 Square-wave approximation of restoring force 
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Fourier transform of displacement 



Fig 2.5 The Fourier transform of the steady state time history 
(Ref: Sundermeyer and Weaver, 1995) 



Fig. 2.6 Spring forces of bilinear and polynomial models 
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CHAPTER 3 


CRACK DETECTION PROCEDURE 
AND NUMERICAL SIMULATION 

The difference between the response of a linear system and that of a nonlinear system, to single 
and multi-tone harmonic excitation is the basis for the proposed crack detection procedure. 
Inherent assumption is that a beam without crack would behave as a linear single-degree-of- 
freedom system, in its fundamental mode, while the presence of a crack would make its stiffness 
bilinear. Further, this bilinear stiffness, approximated through a polynomial, would render the 
system nonlinear. The harmonic components present in the response of a beam, can therefore be 
taken as an indicator of the presence of a crack in the beam. Numerical simulation of such 
systems carried out through MATLAB Simulink models have been discussed in this chapter. 

3.1 Numerical Simulation of the Bilinear Model 

Numerical simulation has been carried to obtain the response of a system governed by the 
following equation. 

mx + cx + (k-k s )x = F 0 smart forx> 0; k e =0 forx< 0 ' (3.1) 

In the above, m,c,k are the mass, damping constant and stiffness of the system, /(/) is the 
excitation force and k e is a stiffness reduction factor introduced to accommodate the crack 
breathing phenomenon. With reference to Fig.2.1, where the crack is on the upper side of the 
beam and the displacement axis, x, is positive downwards, the crack opens for x > 0 and as a 
result there is a reduction in stiffness, denoted by k e * 0 . For the remaining half-cycle, when 
x ^ 0, the crack can be treated as closed and therefore causing no reduction in stiffness 
(*,= 0 ). 

Equation (3.1) can be written in nondimensional form as 

rf + 2 %7j' + (1 - s)rj = sin r r (3 *2) 

where the following parameters are employed 
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natural frequency, p = yfklm ; nondimensional time, x = pt ; 

nondimensional displacement, r/ = x/X st ; damping ratio, £ = c / 2mp (3.3) 

crack severity ratio, s = k s Imp 2 frequency ratio r = a>/ p 

( ) denotes differentiation with respect to r . 

The MATLAB model, employing fourth order Runge-Kutta algorithm, for response simulation 
of equation (3.2) is given in Appendix A. Simulation was initially done for a beam without 
crack, i.e., with e = 0, in equation (3.2) for all 77 . Three sets of excitation 

(i) single tone sinr,r 

(ii) two tone sin r, r + sin r 2 r 

(iii) three tone sin r, x + sin r 2 t + sin r 3 r 

are given to the system. The time and frequency domain plots of these forces are given in 
Figs.3.1 (a)-(c). The response in the three cases is given in Figs. 3.2 (a)-(c). As expected for a 
linear system, the response can be seen to comprise of frequency components at the 
corresponding excitation frequencies, r,,r 2 ,r 3 , alone. 

The response for a bilinear system, with crack severity ratio e = 0.01 , for the three different sets 
of excitation of Figs 3.1 (a)-(c), is shown in Figs. 3.3 (a)-(c). In addition to the frequency 
components at the excitation frequencies, r,,r 2 ,r 3 , the response also shows response harmonics 

at combination frequencies, 2 r, , (r 2 - r, ),(r 3 - 2 r, ), (2 r 2 - r 3 ) and (r, -r 2 +r 3 ) . 

This response pattern is similar to that described in equations (2.24), (2.26) and (2.28), in the 
previous chapter, for a system with polynomial type of stiffness nonlinearity. For, a polynomial 
system, these harmonics are expected to achieve maximum value when the corresponding 
combinational frequency equals the linear natural frequency, p, and since the bilinear system 
exhibits similar response pattern, in Figs. 3.3 (a)-(c), this fact has been employed in the present 
study to explore crack detection. 
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The proposed procedure can be outlined as - 


Step 1 : Identification of the fundamental frequency, p, through a frequency sweep test. 

Step 2: Application of single-tone-excitation, to the system at half of the identified natural 
frequency, i.e. r, = 0.5 . If the system has developed a crack, the response should indicate 
the presence of frequency component 2r,. The excitation frequency, is kept half the 
natural frequency, so as to get, 2r, = 1 , to get the maximum possible magnification for 
this frequency component. This should enable detection of small amounts of nonlinearity, 
which may have developed in the originally linear system. 

Step 3: The presence of nonlinearity can be further confirmed (and eliminating the possibility of 
taking noise as the 2 r, frequency component), by providing a second excitation tone at 
r 2 =1.5 (while keeping r x fixed at 0.5), so as to make (r 2 - r x ) = 1 . This will add a 
magnified ( r 2 - r, ) component to the already existing 2r, component at the natural 
frequency in the frequency response plot. 

Step 4: Further, enhancement of the nonlinear effect can be achieved through addition of a third 
tone of excitation at frequency r 3 = 2.0 . This would add 

(r 3 -2r,),(2r 2 -r 3 )and(r, -r 2 +r 3 ) components to the already existing 2r, and (r 2 -r,) 
components at the natural frequency of the system. 

The numerical simulation based on the above scheme was carried out with the expectation that 
detection of fairly small changes in the system stiffness, consequent to crack development, 
should be possible. The response simulated for four different values of crack severity factors 
e = 0.002, 0.004, 0.006, and 0.008, is shown in Figs. 3.4 - 3.7. The nondimensional excitation 
forces are similar to those described in Fig. 3.1, with frequencies r, = 0.5; r 2 =1.5; r 3 = 2.0 . The 
damping ratio £ is kept constant at 0.001, in all these cases. For each case FFT plots for single 
tone, two tone and three tone excitation have been given. In all cases the presence of a peak at 
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natural frequency (r = 1), with single-tone excitation, indicates the presence of a crack. This peak 
is enhanced for two tone excitation. The enhancement with addition of the third harmonic tone to 
excitation causes only a marginal difference in the peak amplitude at the natural frequency. This 
is explained by the fact that the contribution due to the second and third tone excitation depends 
on the magnitude of the second and third order coefficients, k 2 and k 3 , in the equivalent 

polynomial stiffness representation g(x) = k x x + k 2 x 2 + k 3 x\ It turns out, that for a bilinear 
system, k 3 is a small quantity in comparison to k 2 , resulting in small contributions from the 
third excitation tone (ref. Eqns. 2.22 and 2.28). This is discussed in more detail in the next 
section. It can be readily noted from Figs.3.4-3.7, that as the crack severity ratio e increases, 
crack identification becomes easier due to the presence of larger peaks at the natural frequency. 
The influence of damping ratio on the crack detection process is shown in Figs. 3.8-3.10 and as 
expected for smaller values of damping the response harmonics are magnified making the 
identification process easier. 

3.2 Numerical Simulation of the Equivalent Polynomial Model 

The equivalent polynomial model has been obtained by least square curve fitting of the 
displacement versus spring force plot of the bilinear model. The coefficients k x ,k 2 ,k 3 of the 

equivalent polynomial, resulting from this exercise are employed in the equations (2.24), (2.26) 
and (2.28) to get the response harmonics for single tone, two tone and three tone excitations 
respectively. The expressions obtained from these equations, for the response at natural 
frequency ( r = 1 ), for the three sets of excitation are produced below. 

Single tone ( r x = 0.5 ): 

Z{r = 1) = — H 2 (r, , r, ) + higher order terms 

Two tone (r, = 0.5; r 2 = 1.5 ): 

Z(r = 1) = i H 2 (r, , r, ) + H 2 {~r l ,r 2 )+ higher order terms 
- -k 2 H x (2 r , )[~ H x \r x ) + H { (r 2 )H X (-r, )] 


20 



Three tone ( r, = 0.5; r 2 = 1 .5;r 3 = 2.0 ): 


3 3 

+ - H 3 {r 2 , r 2 ,-r 3 ) + -H 3 (r x -r 2 ,r 3 ) + higher order terms 

T JL 

« -*2 », (>, )[i H- • (r, ) + W, <r 2 )ff , (-r, )] 

3 2£ 2 

+ + r 3 )[—+-{H x (-2r x ) + 2H x (-r x +r 3 )} — * 3 ] 

3 2k 2 

+ -W 1 (r 2 )ff,(r 2 )tf l (r J )/f,(2r 2 - {ff, (2r.) + 2ff,(2r, -r, )}-*,] 

+ 1 H, (r, , (-r 2 )tf , (r, )//, (r, - r 2 + r, ) 

2A: 2 

-r 2 ) + «,(-», +r 2 ) + H,(r, +r, )}-*,] 


(4.6) 


These response harmonics can also be obtained through Runge-Kutta simulation of the 
equivalent polynomial system. The MATLAB model for this purpose is given in Appendix 2. 


It is to be noted that the values of the equivalent coefficients k x ,k 2 ,k 3 are dependent on the range 
of displacement x , in the bilinear model. In the present study, a set of eqivalent polynomial 
coefficients k x ,k 2 , k 3 are obtained for the entire range of displacement x, which are taken to 

remain constant over an entire cycle oscillation. It is possible to divide the displacement range, x 
of the bilinear model, into smaller intervals and do curve fitting for individual intervals, to obtain 
equivalent polynomial coefficients for each displacement interval. This would have increased the 
computational effort significantly. 

The equivalent polynomial model has been created for three crack severity ratios 
(e = 0.002,0.004,0.006). The coefficients k x ,k 2 ,k 3 obtained for each case under single tone 

excitation are 

for e = 0.002 k x = 0.99895081, k 2 = 0.001 18455, k 3 = 0.00014910 

for e = 0.004 k x = 0.99783879, k 2 = 0.00256860, k 3 = 0.00053821 

for * = 0.006 k x = 0.99676018, k 2 = 0.00390345, k 3 = 0.00082413 . 
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Here, damping ratio ^ was taken as 0.001 in all three cases. 

A comparison of the spring forces of the bilinear and equivalent polynomial models is given in 
Figs.3. 1 1 (a)-(c), for the three cases. 

The response of the equivalent polynomial models, for single, two and three tone excitations, for 
the chosen crack severity ratios is given in Figs.3. 12-3. 14. Comparison of these figures with the 
corresponding ones for the bilinear model (Figs. 3.4, 3.5, 3.6) reveals that for lower values of s 
the polynomial model closely represents the response of the bilinear model. For, higher values of 
e, the difference between the two increases. It should be possible to improve upon the 
agreement between the two, through piecewise polynomial approximation, mentioned earlier. 
The present approximation assumes constant polynomial coefficients, through the entire 
displacement range. However, it seems sufficient for a qualitative understanding of the response 
pattern of the bilinear system. 
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(a) single tone excitation 
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(b) Two tone excitation 
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Fig. 3.1 Time and frequency domain plots of excitation forces 









N ondimensional Amplitude Nondimensional AmDlitude Nondimensional Amnlitude 









Nondimensional Frequency 
(b) Two tone excitation (r, = 0.5 ,r 2 = 1.5 ) 



(c) Three tone excitation (r, = 0.5, r 2 = 1.5 ,r 3 = 2.0) 


Fig. 3.3 Response of a beam with crack 
(Crack severity ratio e = 0.01 ) 
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Nondimensionai Amplitude Nondimensional Amplitude Nondimensional Amplitude 





(c) Three tone excitation ( r, = 0.5, r 2 = 1 .5, r 3 = 2.0 ) 

Fig. 3.4 Response of a beam with crack 
(Crack severity ratio e = 0.002) 
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Nondimensional AmDlitude Nondimensional Amplitude Nondimensional Amplitude 



(a) Single tone excitation (r, = 0.5 ) 




(c) Three tone excitation (r, = 0.5,r 2 - 1.5, r 3 - 2 
Fig. 3.5 Response of a beam with crack 
(Crack severity ratio s - 0.004) 





Nondimensional Amplitude Nondimensional Amplitude 




Nondimensional Frequency 
(b) Two tone excitation (r, = 0.5, r 2 =1.5) 



(c) Three tone excitation (r, = 0.5 ,r 2 = 1.5,r 3 = 2.0 ) 


Fig. 3.6 Response of a beam with crack 
(Crack severity ratio e = 0.006) 
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Nondimensional Amplitude Nondimensional Amplitude Nondimensional Amplitude 



Nondimensional Frequency 
(a) Single tone excitation (r, = 0.5 ) 




Fig. 3.7 Response of a beam with crack 
(Crack severity ratio s = 0.008) 
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Nondimensional Amolitude Nondimensional Amplitude Nondimensional Amolitude 





Fig. 3.8 Influence of damping on response of beam with crack (£ - 0.002 ) 
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Nondimensional Amplitude Nondimensional Amplitude Nondimensional Amplitude 



Nondimensional Frequency 
(a) Single tone excitation (r, = 0.5 ) 




Fig. 3.9 Influence of damping on response of beam with crack ( £ - 0.003 ) 
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Influence of damping on response of beam with crack (C = 0.004 ) 
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Fig. 3.1 1 Comparison of spring forces of bilinear and polynomial models 
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Fig. 3.13 Response of equivalent polynomial model ( e = 0.004) 
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Fig. 3.14 Response of equivalent polynomial model ( e = 0.006 ) 
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CHAPTER 4 


EXPERIMENTAL INVESTIGATIONS 

Experimental investigations have been carried out for better understanding of the proposed crack 
detection procedure. The experiments were carried out on a readily available thin beam. Single 
tone, two tone and three tone harmonic excitation signals were generated through Virtual 
Instruments developed in the programming software LabVIEW. Response of the beam was 
obtained for these forces. A transverse crack was then generated in the same beam and the 
exercise was repeated to obtain the response for the three types of excitation. 

4.1 Experimental Set-Up 

The experimental set-up comprised a thin, mild steel beam mounted on an electrodynamic shaker 
at one end (Fig.4.1). The beam approximately represented a cantilever configuration. Two beams 
were chosen for experimentation. The geometric and material properties of the beams are given 
below. 

Beam 1: 

Length, 1=0. 1 70m Width, w=0. 01672m Thickness, d=0. 00054 m 

Modulus of Elasticity, E=200 GP a, Density, p =8000 kg! rr? 

Beam 2: 

Length, 1=0. 170m Width, w=0.01644m Thickness, d=0. 0009m 
Modulus of Elasticity, E=200 GP a, Density, p =8000 kg! n? 

4.2 Instrumentation 

The instrumentation employed for the study is shown in Fig. 4.2. The response of the beam was 
picked up through a miniature accelerometer Type: 4374 Bruel & Kjaer, Charge Sensitivity: 
0.153 pc/ms~ 2 , Voltage Sensitivity: 0.197 mv/ms " 2 . The signal from the accelerometer was 
conditioned and amplified in a Charge Amplifier (Bruel Sc Kjaer, Type 2635). The excitation 
signals for the Electrodynamic Shaker was generated in a Pentium PC, using LabVIEW. 
LabVIEW is program development application, much like C or BASIC. However, LabVIEW 
uses graphical programming language, G, to create programs in block diagram form. In addition 
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of having extensive libraries of functions, LabVIEW includes libraries for data acquisition, GP1B 
and serial instrument control, data analysis, data presentation, and data storage. The programs 
developed in LabVIEW are called virtual instruments (VI s ) because their appearance and 
operation can imitate actual instruments. Thus a VI consists of: 

Front panel : The interactive user interface of the program is called front panel, because it 
simulates the panel of a physical instrument. The front panel can contain knobs, push buttons, 
graphs, and other controls and indicators. Data can be entered using mouse or keyboard. 

Block Diagram: The block diagram is a pictorial solution to a programming problem. It is also 
the source code for the VI 

The Front Panel and Block Diagram for a typical three tone excitation signal, developed during 
the course of this study, are shown in Figs. 4.3 (a), (b). The PC is equipped with an A/D card 
(National Instruments, Texas, AT-MIO-E16, with a scan rate of 250000 samples/s). The signal 
generated by the VI and coming out through the card is taken through an accessory (BNC 2090), 
which connects the A/D card to BNC type connectors, to a Power Amplifier (Type: SS250, MB 
Dynamics, Inc.) The output of the Power amplifier is fed to the electrodynamic shaker. The 
response signal can also be taken to the computer, through the A/D card. Typical VI for data 
acquisition and display is shown in Figs. 4.4 (a), (b). 

4.3 Experimental Results 

Prior to carrying out, crack detection exercise, a Rap-Test was performed on the beam to identify 
its natural frequency. The beam was very slightly perturbed (to keep the vibrations in the linear 
zone) at its free end and the resultant free vibrations were recorded. The FFT of this response is 
shown in Figs. 4.5 (a) and (b), for the two beams. The natural frequency can be seen as 
Beam 1: 14 Hz 

Beam 2: 22 Hz 

The natural frequency were calculated theoretically also, employing the formula 

a>,=W) 2 /WhHP 

where, for a cantilever configuration (/?,/) 2 = 3.52 . The theoretical natural frequencies are 
Beaml: 15.13 Hz 
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Beam 2: 25.21 Hz 

The excitation frequencies, in accordance with the procedure described in the earlier c hap ter, 
were chosen for Beam 1 as: 

Single Tone: co { =7 Hz 

Two Tone: co x = 7 Hz, co 2 = 21 Hz 

Three Tone: ao x = 7 Hz, <a 2 = 2 1 Hz, o> 3 = 28 Hz 

These excitation forces, in the time and frequency domain are shown in Figs.4.6-4.8. 

The FFT of the response, for the beam without crack, for the three sets of excitation are shown 
successively in Figs. 4.9 (a)-(c). It can be noticed that the response harmonics pertain to the 
excitation frequencies, with no combinational harmonics, indicating linearity in the system. 

A transverse notch, of 0.2 mm depth, was then made along the width of the beam, at a distance 
of 127 mm from its free end. This was done by means of a hacksaw, without dismounting the 
beam from the shaker. A Rap Test was again conducted, in order to detect the natural frequency 
of the beam after introduction of crack. No noticeable change in natural frequency was observed 
(from those recorded previously, when there was no crack). The beam was now energised 
through single tone excitation of Fig. 4.6. The response recorded is shown in Fig. 4.10(a). The 
presence of a small 2 u>, = 14 Hz peak, in addition to the one at co x = 7 Hz, corresponding to the 
excitation frequency was an indication of nonlinearity in the system, due to introduction of crack. 
When the two tone excitation of Fig.4.7 was applied, as expected the peak at 14 Hz, was seen to 
increase in magnitude (Fig.4.10 b). The response to three tone excitation of Fig. 4.8 is given in 
Fig. 4.10 (c). No significant change seems to be caused in the peak response at 14 Hz, due to the 
introduction of the third tone in the excitation force. This, as mentioned earlier may be due to the | 
fact that the cubic nature of the nonlinearity may be insignificant in comparison to the square 
term. However, no attempt was made in the present study to determine the nonlinear parameters. 

following set of excitation frequencies were chosen. 

co x = 1 1 Hz j 

co x = 1 1 Hz, co 2 = 33 Hz 1 

co x -\\ Hz, co 2 = 33 Hz, co ^ = 44 Hz 


For Beam 2, the 
Single Tone: 
Two Tone: 
Three Tone: 
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The amplitude for each of the tones was kept identical and equal to those shown in Fig. 4.6-4.8. 
FFT of the response, for the beam without crack, for the three sets of excitation are shown in 
Figs. 4.1 1 (a)-(c). The response after introduction of a transverse notch, of 0.4 mm depth, at a 
distance of 130 mm from the free end is shown in Figs.4.12 (a)-(c). The distinction between the 
response patterns of the uncracked and cracked beam, is not as sharp as in the case of Beam 1 . 
This may be due to an already existing nonlinearity in the beam without crack. Since the 
linear/nonlinear nature is critically dependent on the magnitude of applied force, and the 
consequent displacement level, the need to carry out investigations in a better controlled 
environment was felt. 




Fig. 4.1 Beam with crack mounted on eiectrodynamic shaker 



Fig. 4.2 Instrumentation 
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Fig. 4.3 Description of VI for three tone excitation 
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(a) Front panel 



Fig. 4.4 Description of VI for data acquisition and display 
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Fig. 4.5 Rap test of beams 
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Fig. 4.6 Single tone harmonic excitation for beam 1 



Fig. 4.7 Two tone harmonic excitation for beam 1 



Fig. 4.8 Three tone harmonic excitation for beam 1 
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(c) Three tone excitation(<z;, = !Hz,a> 1 = 2\Hz,oj i = 28 Hz) 


Fig. 4.9 Response of beam 1 without crack 
46 






<D 

I - 

& 

0 -60 



5 20 25 

Frequency (Hz) 

(a) Single tone excitation(co t = 1 Hz) 




Frequency (Hz) 

(b) Two tone excitation^, = 7 Hz,co 2 = 21//z) 



5 20 25 30 

Frequency (Hz) 

(c) Three tone excitation(u), = 7 Hz,q ) 2 = 21 Hz,<d 3 = 28 Hz) 


Fig. 4.10 Response of beam 1 with crack 
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Fig. 4.11 Response of beam 2 without crack 
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Fig. 4.12 Response of beam 2 with crack 
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CHAPTERS 


CONCLUSIONS 


The modeling, numerical simulation and experimentation carried out in this study, suggest that 
multi tone harmonic excitation and higher order FRF analysis can be developed as effective tools 
for transverse crack detection in beams. Bilinear and equivalent polynomial models can account 
for the change in the beam characteristics due to development of a breathing crack, in an 
effective manner. The numerical simulations carried out, also indicate that fairly small bilinearity 
in a system can be detected and magnified through multi tone harmonic excitation. In absence of 
analytical solutions, for a bilinear system, equivalent polynomial models can be built in order to 
rationalise the response patterns. The equivalent polynomial model also helped to outline a step- 
by-step procedure for crack detection. 

The experimental work confirmed the numerically simulated patterns. The present work also 
emphasises the utility of Virtual Instrumentation techniques. A multi-tone excitation force could 
be readily developed in LabVIEW. This would not have been possible through conventional 
means and three electromagnetic shakers with three signal generators would have been required. 

The notch sizes in the experimental beams were, however large, which made crack detection 
possible. The procedure developed can be useful, if small cracks can be detected. There is a need 
to carry out experiments in a more controlled environment. The excitation force needs to be 
carefully controlled. It may be difficult to detect cracks through small excitation forces. Usage of 
higher excitation force levels may however, force beams without crack also, to vibrate in a 
nonlinear regime. Presence of crack may then have to be detected, through comparison of two 
sets of nonlinear response patterns. Experimental work also, needs to be augmented by locating 
the notch at various points along the length of the beam, in addition to varying its size. 

It may be broadly concluded that the technique presented in this thesis, does provide sufficient 
reasons for its further exploration as a tool for crack detection. 
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APPENDIX 1 



SIMULINK MODEL FOR BILINEAR MODEL RESPONSE 



APPENDIX 2 





SIMULINK MODEL FOR EQUIVALENT POLYNOMIAL MODEL RESPONSE 



